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ABSTRACT 

We examine the cosmic star formation rate (SFR) and its dependence on galaxy stellar mass over the redshift 
range 0.8 < z < 2 using data from the Gemini Deep Deep Survey (GDDS). The SFR in the most massive 
galaxies (M* > 10 1() 8 M Q ) was six times higher at z = 2 than it is today. It drops steeply from z = 2, reaching the 
present day value at z~ 1. In contrast, the SFR density of intermediate mass galaxies (10 10 2 < < 10 10 8 M Q ) 
declines more slowly and may peak or plateau at z ~ 1 .5. We use the characteristic growth time tsFR = Pm< / Psfr 
to provide evidence of an associated transition in massive galaxies from a burst to a quiescent star formation 
mode at z ~ 2. Intermediate mass systems transit from burst to quiescent mode at z ~ 1, while the lowest mass 
objects undergo bursts throughout our redshift range. Our results show unambiguously that the formation era 
for galaxies was extended and proceeded from high to low mass systems. The most massive galaxies formed 
most of their stars in the first ~ 3 Gyr of cosmic history. Intermediate mass objects continued to form their 
dominant stellar mass for an additional ~ 2 Gyr, while the lowest mass systems have been forming over the 
whole cosmic epoch spanned by the GDDS. This view of galaxy formation clearly supports 'downsizing' in 
the SFR where the most massive galaxies form first and galaxy formation proceeds from larger to smaller mass 
scales. 

Subject headings: Galaxies: formation — Galaxies: evolution 



1. INTRODUCTION 

The evolution of the global star formation rate provides a 
sensitive probe of galaxy formation and evolution. The ear- 
liest determinations of the evolving star formation rate den- 
sity (SFRD) showed a steep decline from z ~ 1 to the present 
(Lilly et al. 1996; Madau et al. 1996). The behavior of the 
SFRD at early epochs (z > 1) remains uncertain due to vari- 
ations amongst SFR diagnostics and poorly constrained (yet 
potentially large) extinction corrections in the primary rest- 
frame UV diagnostics (Steidel et al. 1999). Despite these 
challenges Hopkins (2004) recently compiled results from 33 
studies over the range z = to z ~ 6. The data were used 
to constrain the luminosity function of star-forming galaxies, 
and were found to be consistent to within a factor of three over 
0<z<6. 

Additional insight into star formation histories (SFH) may 
be gained by the consideration of other physical galaxy prop- 
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erties. Stellar mass is arguably the key parameter. The mass 
in stars provides a measure of the integral of past galaxy 
stellar mass assembly, which can be coupled with the in- 
stantaneous SFR to give a more complete view of galaxy 
evolution. Mass-based evolution studies are far more de- 
terministic since unlike luminosity, mass evolution is mono- 
tonic. While high-luminosity galaxies often evolve into low- 
luminosity systems, massive galaxies at early epochs must 
have descendants among the present massive galaxy popula- 
tion. Recent advances in the modeling of multi-color, and 
particularly near-IR, selected samples lead to fairly robust de- 
terminations of stellar masses for galaxies over a wide range 
of redshifts and luminosities (Brinchmann and Ellis 2000; 
Fontana et al. 2004; Glazebrook et al. 2004, hereafter Pa- 
per III). 

Heavens et al. (2004) inferred the star formation history 
(SFH) of the universe by modeling the spectra of 96,545 local 
SDSS galaxies. Their results indicate that SFHs vary strongly 
with present-day stellar mass. Galaxies with M* < 10 10 ' 7 M Q 
go through the peak of their star formation (SF) at low red- 
shift (z < 0.5) and > 10 11 2 M Q galaxies show an increas- 
ing SF activity to z > 2. We take advantage of the near-IR 
selection of the Gemini Deep Deep Survey (GDDS, Abraham 
et al. 2004, hereafter Paper I) to study directly the evolution 
of the SFRD as a function of stellar mass at the epoch of ob- 
servation. Direct measurements of SFRs in the most massive 
systems at z > 1 were made possible only recently, with the 
progress of spectroscopy of near-IR surveys probing the red- 
shift range 1 < z < 2. We use (Hq, VLm, ^a) = (70 km s" 1 
Mpc" 1 , 0.3, 0.7) and Vega magnitudes throughout this paper. 

2. DATA AND SAMPLE SELECTION 

The GDDS is a spectroscopic survey of an optical and near- 
IR selected sample targetting massive galaxies at 0.8 < z < 
2.0. The sample is drawn from the Las Campanas IR imag- 
ing survey (McCarthy et al. 2001; Chen et al. 2002), and 
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was designed to select galaxies in the 0.8 < z < 2 range with 
an emphasis on the reddest galaxies. Galaxy stellar masses 
are determined from the mass-to-light ratio M/Lk obtained 
by fitting the VIK photometry with a grid of SED models (Pa- 
per III). The sample selection function and weights, details of 
the observations and catalogs are described in Paper I. The 
sampling weights are derived as a function of color and mag- 
nitude and are used when computing volume-averaged quan- 
tities such as the mass density (Paper III) and the SFRD. 

The sample selected to compute the SFRD consists of spec- 
tra from the GDDS that satisfy: (1) K < 20.6 (survey limit), 
(2) redshift confidence level greater than 75% (conf > 2 in 
Paper I notation) and (3) absence of strong active galactic nu- 
clei activity (agn = in Paper I notation). From the original 
GDDS sample of 308 spectra, 211 meet the ^-detection and 
redshift confidence class criteria. Strong AGN contamination 
occurs for 1.9% (4/211) of those objects, bringing the final 
sample to 207 galaxies. 

3. METHOD 

Given the redshift range spanned by the GDDS, the avail- 
able SFR indicators are the [OII]]A3727 emission line and the 
luminosity of the rest-frame UV continuum. For the latter, we 
chose the absolute rest-frame AB magnitude M2000, defined 
in a synthetic 1900 < A < 2100 A box filter using an em- 
pirical interpolation scheme from the observed V and / mag- 
nitudes (Savaglio et al. 2004 (Paper II)). The redshift range 
is restricted to z < 1.6 for SFR([OII]) measurements as the 
[Oil] emission is redshifted out of the optical range at z > 1 .6. 
For the SFR(M2ooo) measurements, we use the redshift range 
where the M2000 interpolation is reliable, i.e. 1.2 < z < 2.0. 
Details of the SFR measurements will be given in a forth- 
coming paper (Juneau et al. in prep.). Briefly, we use the 
standard SFR(//a) conversion of Kennicutt (1998) assuming 
([Oll]/Ua) obs = 0.5 (Glazebrook et al. 1999). Since we use 
the observed ratio, we need to correct the integrated luminosi- 
ties at Ha. We adopt the average extinction of An a = 1 derived 
from local galaxy samples (Kennicutt 1992). 

For SFR measurements based on the rest-frame UV contin- 
uum, we apply a dust attenuation correction of the stellar con- 
tinuum of Ay = 1 • This is the value typically used in the liter- 
ature when no direct measure of dust obscuration is available 
(e.g. Lilly et al. 2003, Sullivan et al. 2000). Following the 
prescription of Calzetti (2001) this correction is A2000 = 2.2 
at 2000 A. The mean dust obscuration in galaxies depends 
on sample selection, /^-band selection could include more 
heavily obscured systems. If the mean extinction is greater by 
1 mag, the values of SFR (§4) and SFRD (§5) will shift up 
by a factor 2.5 whereas tsFR (§5) will be lower by the same 
factor. 

The SFRD is computed with the 1 /V max method and cor- 
rected for both sampling and spectroscopic incompleteness. 
We define spectroscopic completeness factors that depend on 
the color and the magnitude in the same fashion as the sam- 
pling weights from paper I. We divide the number of spec- 
tra with high confidence redshifts (conf > 2) by the number 
of slits in each cell of the color-magnitude plane (see Figure 
12 of Paper I). The median spectroscopic completeness of the 
sample according to this definition is 82%, consistent with the 
overall spectroscopic completeness of 79% for the GDDS. 

Throughout this paper, we adopt the Baldry & Glazebrook 
(2003) (BG03) Initial Mass Function (IMF), which has a very 
similar slope to the Salpeter (1955) IMF at high masses and 
provides a good fit to cosmic and galaxy colors locally. The 



galaxy masses and SFRs based on BG03 can be used in- 
terchangeably with those of Salpeter given the conversion 
M(BG03) = 0.55M(SP), this ratio is virtually independent of 
SFH to an accuracy of a few percent. 

4. STAR FORMATION RATES 

The SFRs of the individual galaxies are plotted in Figure 1 . 
The values estimated from the [Oil] luminosity (circles) and 
those obtained with the continuum luminosity at 2000A (tri- 
angles) suggest an increase by over one order of magnitude 
in the upper envelope of the SFR values between z = 0.8 and 
z = 2. The color of the plotting symbols is keyed to the stel- 
lar mass of the galaxies. The mass bin corresponding to the 
lowest mass galaxies is 10 9 < < 10 10 - 2 M© (blue), the in- 
termediate mass is defined as 10 10 2 <M t < 10 10 ' 8 M Q (green) 
and the high mass galaxies have 10 10 ' 8 < M* < 10 1L5 M Q 
(red). Note that /T-selected samples will miss contributions 
from low-mass star-forming galaxies fainter than the Zf-limit. 
(We show later in Figure 2 that our full K < 20.6 sample un- 
derestimates the total SFRD by a factor of 2-3.) 
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FIG. 1.— SFR derived from L([OH]) (circles) and from L(2000A) (tri- 
angles). The values are corrected for obscuration by dust using An a = 1 
for emission-line measurements and Ay = 1 for UV continuum measure- 
ments (A2000 = 2.2). The open symbols show objects with K s > 20.6 and 
the symbols color is keyed to the mass bin in which the object belongs: 
1O 9 M < M, < 10 10 2 Mq (blue); 10 la2 M Q <M„ < 1O 1O ' 8 M (green) and 
1O 1O 8 M <M* < 1O 1L5 M (red). The cases where [Oil] emission is not 
formally detected are illustrated with upper limit symbols. The objects with 
no formal [Oil] detection allow us to estimate our [Oil] flux sensivity limit 
< 0.6 X 10~ 17 ergs~'cirT 2 (dotted line). A subsample of X-ray star-forming 
galaxies from HDF-N and Chandra data are shown as a comparison (plus 
symbols, from Cohen (2003)). The points outside our target redshift range 
(gray area) are not used in subsequent analysis. 

In nearly 10% (20/207) of the galaxies, no [Oil] emission 
is detected in our spectra. These systems appear to be primar- 
ily massive, quiescent galaxies and their number is uniformly 
distributed with z (upper limit symbols). The [Oil] flux limit 
varies from one spectrum to another, depending on the red- 
shift and the integration time for each mask. A conservative 
flux detection limit of F([OH])ij m = 0.6 x 10" 17 erg s _I cm~ 2 is 
shown as a dotted line in Figure 1 . The upper envelope of the 
[Oil] -derived SFRs is representative of star formation in mas- 
sive galaxies selected at longer wavelength (e.g. 2/im) rather 
than samples selected in the rest-frame UV. To contrast our 
K selection, Figure 1 also shows a complementary dataset: a 
subsample of HDF-N galaxies detected by Chandra (shown as 
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plus symbols, and taken from Cohen 2003). Since strong X- 
ray and [Oil] emission is necessary for inclusion in this sam- 
ple it includes objects with very high SFRs and the sample is 
not directly comparable. However, it follows the same trend 
as the GDDS sample. Figure 1 demonstrates that the photo- 
metric redshift cut used to select galaxies at 0.8 < z < 2.0 was 
efficient, as only a handful of objects were observed outside 
that redshift range (gray area). 

The two tracers of SFR were directly compared in the red- 
shift interval in which they overlap: 1.2 < z < 1.6. We found 
a good linear correlation, with a dispersion ~ 0.6 dex (Juneau 
et al. in prep.). The scatter might result from the fact that 
adopting a single value for dust attenuation is not likely to be 
a realistic representation of a diverse set of galaxies. 

5. STAR FORMATION RATE DENSITY AND CHARACTERISTIC 
GROWTH TIMESCALE 

The mass selection of the GDDS provides an opportunity 
to probe cosmic star formation as a function of galaxy stel- 
lar mass. In Figure 2 we show the global SFRD for our three 
mass bins, the total GDDS sample, and the compilation from 
Hopkins (2004). As in Figure 1, the symbol type is keyed to 
star formation indicator, circles show [Oil] measurements, the 
triangles represent rest-frame 2000A based rates. The local 
SFRDs derived from the SDSS by Brinchmann et al. (2004) 
in our three mass bins are also shown in Figure 2 (square sym- 
bols). These have been transformed to the same IMF and cor- 
rected for extinction in the same manner as the GDDS mea- 
surements. The comparison with Figure 1 indicates that the 
increase of SFRD with z is linked with the increase of SFR in 
individual galaxies. 

Paper III shows that for > 10 10 8 M Q our sample is 
mass-complete over our full redshift range. At lower masses, 
10 10 - 2 < < 1O 1OS M we start to become incomplete for 
z > 1.15. However the effect is minor because it is the red ob- 
jects which start to be missed and these contribute least to the 
SFRD. In this sense, mass-incompleteness also means sensi- 
tivity biased towards bluer objects. We can estimate this in- 
completeness, and calculate a reasonable correction, by boot- 
strapping from the z < 1.15 complete sample. We calculate 
that if these objects were placed at z = 1 .5 we would see 90% 
of the total SFRD, and at z = 1 .8 we would see 62%. Note we 
are effectively assuming that the high-z sample displays the 
same range of M/L K values as the low-z sample. In practice 
we expect the high-z sample to be bluer, which would make 
us more complete. Thus we think our correction is conser- 
vative; we extend the upper error bars in Figure 2 to reflect 
the magnitude of this correction. For M* < 10 1() 2 M© we are 
incomplete throughout our redshift range and plot the SFRD 
values as lower limits. 

The key result of this paper is illustrated in Figure 2: the 
cosmic SFRD is a strong function of galaxian stellar mass. 
The SFRD in our high mass bin (M t > 10 10 8 M Q ), while mak- 
ing a minority contribution to the global SFRD, is a factor ~ 6 
higher at z = 2 than it is at present. The SFRD of these mas- 
sive galaxies has strongly declined since z = 2 and reached the 
present-day level at z ~ 1 . The SFRD in the intermediate mass 
bin (10 10 - 2 < M* < 10 10 8 M Q ) has also steeply declined since 
Z ~ 1 .2 and appears to have reached a peak or a plateau at a 
redshift of ~ 1 .5. The SFRD in the intermediate mass bin then 
declined to the present day level at z < 1 and since that time 
most of the SFRD has been in low mass galaxies. Our lowest 
mass bin (M* < 10 10 2 Mq) is strongly mass-incomplete. The 
comparison with the Hopkins points shows that we miss about 



half the total SFRD at z ~ 1 . This would come from low mass 
galaxies fainter than our survey K limit. 
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FIG. 2.— SFRD derived from L([OII]) (circles) and from L(2000A) (tri- 
angles). The symbols are color-coded with the mass ranges as in Figure 1. 
The error bars in redshift show the width of the redshift bins used. The er- 
ror bars in SFRD combine shot noise and mass-completeness corrections. 
Both the sampling and the spectroscopic completeness corrections were ap- 
plied. The squares are the values found locally by Brinchmann et al. (2004) 
converted according to our assumed IMF and dust correction. The com- 
pilation made by Hopkins (2004), where all the values are converted to a 
(Q.M = 0.3, f2y\ = 0.7, ft = 0.7) cosmology, are overplotted with diamonds. The 
line is the fit derived by Cole et al.(2001) assumming Ay = 0.6 (solid line). 

Additional insight into the growth of galaxies can be gained 
from comparing the stellar mass density to the SFRD over 
a range of redshifts. The ratio gives a characteristic growth 
timescale tsFR = Pm, I Psfr which can be interpreted simply as 
the time required for the galaxies to assemble their observed 
stellar mass assuming that their observed SFR stays constant. 
Results for the GDDS sample are shown in Figure 3. In this 
figure we compare tsFR with the Hubble time f#(z), the age of 
the universe at a given redshift. If equal, galaxies can form 
all their observed stars in a Hubble time. At a given redshift, 
tsFR > tH suggests that the galaxies are in a declining or quies- 
cent star formation mode at the observed redshift and that the 
bulk their star formation has occurred in the past. Conversely, 
hFR < tH indicates that the galaxies are going through a burst 
phase. This past average-to-present SFR allows us to investi- 
gate the mode of SFR in the SFRD(z) diagram. The highest 
mass galaxies are in a quiescent mode at low-z, transiting to 
a burst mode at z ~ 1.8. The intermediate mass objects make 
a similar transition at z ~ 1.1. The lowest mass systems ap- 
pear to be observed in a burst mode at all redshifts, though 
we note our selection is only sensitive to the bluest objects for 
this mass-incomplete subsample. The transition redshift from 
burst mode to quiescent mode is a strong function of stellar 
mass and looks like a time-delayed echo of the corresponding 
downturn in the SFRD(z) diagram. 
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FIG. 3. — Characteristic timescale of stellar mass growth in galaxies. The 
SFRD is derived from L([OII]) (circles) and from L(2000A) (triangles). The 
plotting symbols are keyed to the galaxy stellar mass as in Figures 1 & 2. 
The error bars in redshift show the width of the redshift bins used. The error 
bars in SFRD are statistical. The gray line shows the age of the universe in 
our adopted cosmology. It indicates the transition from quiescent SF mode 
to burst SF mode (gray area). Along the dotted lines, lookback time to the 
present allows galaxy stellar mass to increase by a factor of 2, 5 or 10, if the 
SFR stays constant, as labeled. 



rising SFRD seen in the high-mass bin of Figure 2 continues 
to increase, reaching a peak somewhere in the 2 < z < 4 inter- 
val. The most massive local galaxies seen in the SDSS also 
appear to be dominated by stars with early formation redshifts 
(Heavens et al. 2004). The Hubble sequence begins to emerge 
between 1 < z < 2, and appears to be in place at z ~ 1 (Brinch- 
mann & Ellis 2000). Disk galaxies that are seen at z > 2 are 
primarily in massive stellar systems (e.g. Labbe et al 2003; 
Stockton et al. 2004). As Figure 3 shows, by z ~ 1 both the 
high and intermediate mass populations, those that dominated 
the elliptical and disk portions of the Hubble sequence, have 
transitioned to fairly quiescent SF and have formed the bulk 
of their stars. A puzzle remains in the continued accumula- 
tion in the number of massive (see Paper III and Fontana et 
al. 2004) and/or red-sequence (Bell et al. 2004) galaxies by a 
factor of ~ 2 over < z < 1 . This now has to occur without 
making a large contribution to the global SFRD in the visible 
and rest-UV, perhaps via mergers that produce either little or 
heavily obscured star formation. 

Our view of galaxy formation reveals a clear 'downsizing' 
of star formation from high-mass to low-mass galaxies with 
time. This basic picture was first suggested by Cowie et al. 
(1996). The downsizing picture of galaxy stellar mass assem- 
bly is the reverse of the original picture of hierarchical dark 
matter assembly (Blumenthal et al. 1984) where the large 
units come last in the sequence. If the hierarchical picture is 
correct, then star formation must be much more efficient in 
early times in high-mass systems, as is required (e.g. Paper 
III) to explain their space densities at high redshift. 



6. DISCUSSION AND SUMMARY 

Figures 2 and 3 paint a simple picture of galaxy formation 
in a volume-averaged sense. The upper envelope of SFR rises 
with redshift and one sees very clear mass-dependent effects 
in the SFRD(z) and tsFuiz) diagrams. The mass-scale of the 
SFR and of burst activity both decline with time. 

The basic conclusion is that the most massive systems 
formed early and were nearly finished forming their stars by 
z ~ 1.5-2. This is supported by several other lines of evi- 
dence. The reddest galaxies in the 1 < z < 2 interval appear to 
have ages > 1 Gyr and z/ > 2 (McCarthy et al. 2004; Cimatti 
et al. 2004). These objects are ~ 10 n M Q and suggest that the 
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